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Abstract 

The space of local operators in three-dimensional quantum electrodynamics contains 
monopole operators that create n units of gauge flux emanating from the insertion point. 
This paper uses the state-operator correspondence to calculate the anomalous dimensions 
of these monopole operators perturbatively to next-to-leading order in the 1/Nf expansion, 
thus improving on the existing leading order results in the literature. Here, Nf is the number 
of two-component complex fermion flavors. The scaling dimension of the n — 1 monopole 
operator is 0.265A/ — 0.0383 + 0(1 /Nf) at the infrared conformal fixed point. 
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1 Introduction 



There are several reasons why one may be interested in studying quantum electrodynamics 
in three dimensions (QED 3 ). To a high energy physicist, QED 3 may be of interest be- 
cause it bears some resemblance to four- dimensional quantum chromodynamics (QCD 4 ). 
Indeed, both theories are asymptotically free, and they are both believed to exhibit spon- 
taneous "chiral" symmetry[j] breaking when the number of fermion flavors in QED3 is suffi- 
ciently small [l}j3]. At the same time, QED 3 is in many ways simpler than QCD 4 . While 
QCD 4 is a non-abelian gauge theory that requires non-trivial renormalization and in which 
asymptotic freedom follows from a complicated beta-function computation [4^5], QED3 is a 
super-renormalizable abelian gauge theory where asymptotic freedom follows from dimen- 
sional analysis. In three dimensions, the gauge coupling e 2 has units of mass, and hence 
the dimensionless coupling e 2 //x, /i being the renormalization group (RG) scale, becomes 
arbitrarily small in the ultraviolet. Therefore, understanding the dynamics of QED3 should 
be in principle less onerous than understanding that of QCD4, and at the same time QED3 
may provide us with valuable lessons about the physics that leads to spontaneous chiral 
symmetry breaking. 

To be precise, let Nf be the number of two-component complex fermions in QED3. The 
following discussion will be restricted to the case where Nf is even, because otherwise it would 
be impossible to regularize the theory while preserving parity and time reversal symmetry. 
In Euclidean signature, the action is 



S = d 3 r 



N f 

-i^i^ + ^t (ZIP + 4) 



4e 2 

a=l 



'1.1: 



where ip a are the fermion fields, is a U(l) gauge field with field strength F^ u , and Ip is 
the Dirac operator. This theory has a global U(Nf) symmetry under which the fermions 
ip a transform as a fundamental vector. Chiral symmetry breaking refers to the breaking of 
this symmetry to U (Nf/2) x U (Nf/2). While the maximal value N^ nt for which one expects 
this symmetry breaking to occur has been the subject of some debate, recent lattice gauge 
theory results |6]j8j suggest that Nf nt = 2 J 2 ] For Nf > 2 the theory is believed to flow to a 



lr rhere is no chirality for fermions in three dimensions. The symmetry in question would be a chiral 
symmetry if the same theory were considered in four dimensions. 

2 Other approaches give higher critical values Nj Tlt . For instance, an analytical approach based on solving 
the Schwinger-Dyson equations self-consistently gives N c f rit = 6 m\n\. Bounds obtained by using the F- 

'?" ' ' ' 



theorem 12f|15 give N c f rit < 6 16 
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strongly interacting conformal field theory (CFT) in the infrared. At least at large enough 
Nf, the CFT is obtained by simply taking the limit e 2 — > oo in the action (1.1). 

To a condensed matter physicist, QED3 may be of interest because it describes, for 
instance, the effective low-energy dynamics of "algebraic spin liquids" jl7|. As explained 
in fl8], SU(Nf) spins on a lattice can be described in a slave fermion formalism 19 as 



compact QED3 with Nf fermion flavors. Compactness of the U(l) gauge field means that 
monopole configurations are not suppressed. Just like in the case of no fermions studied 
by Polyakov 20 21 , at small Nf one can argue that these monopoles proliferate, and their 
proliferation leads to confinement of electric charges. In renormalization group language, 
the monopoles can proliferate provided that the operators that create them are relevant 



in the RG sense 22,23 . When Nf is large, it can be shown that the monopole operators 
become irrelevant 24 , and therefore one can ignore the compactness of the gauge group 1 22 1 . 
All other relevant operators, such as fermion bilinears and Chern-Simons interactions, are 
suppressed provided that one requires invariance under various discrete symmetries, in which 
case the resulting low-energy effective theory, namely the CFT obtained by taking e 2 — > 00 



in (1.1), is a stable RG fixed point. This stable CFT is the low-energy effective theory of an 
algebraic spin liquid. 

In this paper, I will calculate the scaling dimension of the monopole operator that inserts 
one unit of magnetic flux by performing an expansion to next-to-leading order in 1/Nf. 
The leading term in this expansion was found by Borokhov, Kapustin, and Wu [24]. I 
will therefore improve on their result. While the validity of the 1/Nf expansion has not 
been tested in QED3, there is evidence that in supersymmetric versions of QED3 the 1/Nf 
expansion gives reasonably accurate results even for small values of Nf 25,26 . The tests 
of the 1/Nf expansion were made possible in supersymmetric theories by the technique of 
supersymmetric localization 27 28 combined with F-maximization 29,30 , whereby one can 
obtain exact results for various supersymmetry-protected quantities, even at strong coupling. 
It is therefore desirable to make use of the 1/Nf expansion in non-supersymmetric theories 
as well, and this paper will focus on QED3. 

In QED 3 , the scaling dimensions of the operators in the zero- monopole sector can be 
computed in the 1/Nf expansion using Feynman diagrams, and many results are known to 
several orders in 1/Nf (see, for example, [TtJ ) . The situation is significantly harder when one 
includes monopoles, because the monopole creation operators do not have simple expressions 



in terms of the fundamental fields of the theory (1.1), so conventional Feynman diagrams 
are not of much use. The monopole scaling dimensions can be computed, however, using the 
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state-operator correspondence J^] which maps local operators of a CFT inserted at the origin 
of IR 3 to normalizable states on S 2 x M, the IR coordinate being interpreted as Euclidean 
time. The scaling dimension on IR 3 is mapped to the energy of the state on S 2 . 

A monopole operator of strength n is a local operator that changes the boundary condition 
for the gauge field at the insertion point such that the field strength F = dA integrates to 

F = 2im (1.2) 

s 2 

over any (sufficiently small) two-sphere surrounding the insertion point. Here, Dirac quanti- 
zation imposes n G Z. The corresponding state on S 2 is the ground state in the presence of 
n units of magnetic flux through S 2 ; the excited states in the presence of this magnetic flux 
correspond to composite operators that contain a monopole operator of strength n. Standard 
thermodynamics equates the ground state energy with the free energy on S 2 , which can be 
computed as minus the logarithm of the S 2 x M. partition function. The scaling dimension 
of a monopole operator of strength n can therefore be extracted from the S 2 xK partition 
function in the background of n units of magnetic flux through the S* 2 . 

The partition function on S 2 x IR can be evaluated perturbatively in 1/Nf as follows. 
To leading order in Nf one can ignore the fluctuations of the gauge field and perform a 
Gaussian integral over the fermionic fields. The results obtained by this procedure agree 



with those of 24 . The next order term in the large Nf expansion comes from performing 
a functional integral over the fluctuations of the gauge field. This integral is also Gaussian, 
because the higher order terms in the effective action for the gauge field fluctuations are 
suppressed by positive powers of 1/Nf. The evaluation of this Gaussian integral is the main 



focus of this paper. Similar but significantly simpler computations were performed in 33 
where the contribution to the S 2 ground state energy coming from a scalar fluctuation was 



computed in a similar setup, and in 25 , where an integral over gauge field fluctuations on 
S 3 was computed in the absence of any monopoles. 

In contrast with the computation presented in this paper, a similar computation in a 
supersymmetric theory would be much simpler if, as mentioned above, one appeals to super- 



symmetric localization and F-maximization 29 . Indeed, in M = 2 supersymmetric QED 



with N chiral super-fields with charge +1 and N chiral super-fields with charge —1 under the 
£7(1) gauge group, the scaling dimension of a BPS monopole operator with strength n = 1 



3 For other instances where the scaling dimensions of monopole operators were computed using the state- 
operator correspondence, see 24pT ■ 33 . See also 34 for a different method of computing scaling dimensions 
of monopole operators. 
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is N(l — A) (see for example [26]), where A is the scaling dimension (or R-charge) of one of 



the chiral super-fields. This scaling dimension was computed using F-maximization in 25 



both exactly and in a 1/N expansion, and the two calculations agree very well even at small 



values of iV — see Figure 2 in 25 . The scaling dimension of the n = 1 monopole operator is 
in this case 

N 2 2(tt 2 - 12) ^ /Ar 2x n q\ 
1 1 - + 0(N~ ). (1-3) 

2 7T 2 7l 4 N { ' V ' 

In the non-supersymmetric case presented below, there are no methods that yield a simple 



analytic answer like (1.3), and one has to resort instead to numerical methods. 

The remainder of this paper is organized as follows. Section [2] contains the setup of 
the problem in more detail. In Section [3j I reproduce the leading order results of 24 



by performing the Gaussian integral over the fermionic fields. In Section |4j I outline the 
strategy used for computing the leading l/Nf correction to these results. Section [5] contains 
expressions for the fermion Green's function on S 2 x R in the presence of n units of magnetic 
flux through S 2 . This Green's function is the central ingredient in the effective action for 
the gauge field fluctuations. In Section |6j I evaluate the Gaussian integral over the gauge 
field fluctuations. Lastly, Section [7] contains a discussion of these results. 

2 Setup and conventions 

On an arbitrary conformally flat Riemannian three-manifold, the QED 3 action with Nf 
two-component complex fermions is 

(2.1) 



f 

£1 = 1 J 



where g is the determinant of the metric, ip a are the Nf fermion fields, and the gauge field is 
written as the sum of a background and a small fluctuation around this background. 
We are interested in studying this theory on S* 2 x K in the background of n units of magnetic 
flux through S 2 . Parameterizing S 2 x R using coordinates r = (8, <ft, r) such that the metric 
is written as 

ds 2 = dd 2 + sin 2 Odcf) 2 + dr 2 , (2.2) 
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one can take the background gauge field to be 



A{r) = -(1 -cos 9)d(j). (2.3) 



The field strength J 7 = dA integrates to 2im over S 2 , as in (1.2). The expression (2.3) is 
well-defined everywhere away from the South pole at 6 = ir, where there is a Dirac string 
extended in the IR direction. The requirement that this Dirac string should be invisible 
restricts n G Z. We will assume n > 0. 

In working with spinors on a curved manifold, one should specify the conventions used 
for the frame and gamma matrices. To simplify the subsequent analysis, we can introduce 
the frame obtained from conformal transformation of the standard one in IR 3 : 

e l = e~ T dx\ x = e T x = e T (sin 6 cos <p sin # sin cos6>) , (2.4) 

which can be found by writing the standard line element on IR 3 in spherical coordinates as 



dx 2 = dp 2 + p 2 (d6 2 + sin 2 ddcf 2 ) such that the metric on S 2 x IR is 



ds 2 = %, p = e. (2.5) 
p 2 

We will use the gamma matrices 7, = Oi where Oi are the Pauli matrices]^] 

The following sections are devoted to the computation of the ground state energy = 
- log on S 2 x IR in the presence of the background (|2.3l). This quantity can be expanded 
at large Nf as 

FW = NfF™ + Fi n) + .... (2.6) 

One expects F^ = because when n = the ground state on S 2 corresponds to the identity 
operator on IR 3 , which has vanishing scaling dimension. While we will check explicitly that 
Fq = 0, we will take F^ = as an assumption and identify F^ with F^ — F±°\ 



3 Leading order free energy 

To leading order in Nf we can ignore the fluctuations of the gauge field, and the action 



becomes that of free fermions in the background gauge field (2.3). This action can be 



4 There is no difference between upper and lower frame indices in Euclidean signature. 
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written in almost diagonal form by expanding the fermion fields in terms of Fourier modes 
in the ffi. direction as well as in terms of analogs of the S 2 spherical harmonics that are 



appropriate for describing a spin-1/2 charged particle in the monopole background (2.3) 



The main ingredients in constructing these harmonics are the monopole spherical harmonics 



Y, 



n/2,£m 



35 36 , with £ > 1/2 and — t < m < £, which are simultaneous eigenfunctions of the 



gauge- covariant angular momentum operators: 



L 2 Y r 



n/2,tm 



+ l)Y n/ 



21m 



L»Y r 



2,ln/2,im 



mY r 



n/2,£m 



(3.1) 



For explicit formulas for Y n / 2 / m m terms of the angles on S 2 , see 35 36 or Appendix [Aj The 
monopole spherical harmonics form a complete basis of functions on S 2 suited for describing 



spinless particles with unit gauge charge in the background of the monopole (2.3). When 
n = they are nothing but the usual S 2 spherical harmonics. 



As explained in more detail in Appendix A, the scalar harmonics Y n n^ m can be general- 
ized to include spin by simultaneously diagonalizing S 2 , L 2 , J 2 , and J3, where S, L, and J 
are the spin, orbital, and total angular momentum operators. The spin-1/2 case yields two 
sets of spinors that we denote by S n / m and T n ^ m . They have orbital angular momentum 
quantum number i and total angular momentum quantum numbers 



+ 



l- 



mi 



m 



m + 



(3.2) 



The expansion of the fermion fields reads 



2^ 



E ^T(^)Tn,im(0, 



E 



, (3.3) 



where and I'Jj are (anti-commuting) coefficients. In this expression and throughout 
this section all spinor indices carried by i/; a , S n< £ m , and T n ^ m are suppressed. The range of 
m in the sums in (3.3) follows from (3.2) and the usual —j < rrij < j. While the spin-1/2 
monopole harmonics S n ^ m and T n ^ m diagonalize S 2 , L 2 , J 2 , and J 3 , they are not eigenspinors 
of the gauge-covariant Dirac operator ilp + Jl because this operator does not commute with 
L 2 . Indeed, for given j = I — 1/2 and rrij = m + 1/2, there are only two spinors that differ 
in their L 2 quantum number, namely T n ^_i^ m and S n ^ m ; starting from 24 it can be shown 
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that 



^n % £m& 



N n , (cj + iM r 



n,(t-~L)m e 



(3.4) 



where the matrices M n ^ and N n ^ are given b)j^] 



M 



n,£ 



4£ 2 



4£ 2 



4<? 2 



1 4£ 2 



N 



21 



1 4£ 2 



1 - 21 



u 

2( 



(3.5) 



Note that these matrices act trivially on the spinor indices of T n u_ \\ m and S n ^ m , which, as 
mentioned above, are consistently being suppressed. 

The case i = n/2 deserves a comment. Since there are no scalar monopole spherical 
harmonics with orbital angular momentum less than n/2, when I = n/2 (or equivalently 



when j = (n— l)/2) the matrices in (3.5) should be thought of as 1 x 1 matrices equal to the 
bottom right entries of the expressions in (3.5). In particular, M„ in / 2 = and N ni „,/ 2 = 1. 



The modes S n i n /2)m are the n zero modes of the Dirac operator on S 2 in the presence of n 
units of magnetic flux. 



Using (3.4) and (3.3) and ignoring the gauge field fluctuations, the action takes the 



block-diagonal form 



N f r a 00 l ~ l At/( £ - 1 )"V ■ 

a=X J £=^m=-t \ ^a,S 



CO 



■ (3-6) 



After performing the Gaussian integral over the fermions, the logarithm of the S 2 x 
partition function is thus 



/, oo 
— ^2llogdet (N q/ (u + iM q>t )) 
n e= q 



(3.7) 



the factor of 2t coming from the sum over m, or equivalently from the 2j + l = 21 degeneracy. 
The determinant in (3.7) is easily computed, and the coefficient in (2.6) can be expressed 
as 



u 2 + f 



n 



4 i (3-« 
5 I thank Mark Mezei for helping me correct a sign error in a previous version of this equation. 
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n 


^0 








1 


0.265 


2 


0.673 


3 


1.186 


4 


1.786 


5 


2.462 


6 


3.206 



Table 1: The values of Fq as obtained in a few particular cases by evaluating (3.11). 



This expression is divergent and requires regularization. One way of extracting the finite part 
is to write log A = —dA~ s /ds , followed by an evaluation of the sum and integral in (3.8) 



at values of s where they are absolutely convergent, and then by an analytic continuation of 
the answer to s = 0. (See also |33|, where a similar expression was regularized in the same 
fashion.) Performing the u integral yields 



0— n 
c ~ 2 



n 

T 



(3.9) 



After adding and subtracting quantities that are divergent when s = 0, one can write this 
sum as 



- 2 E 



n 



p-2 S + !L(i_ 2s )r 2s 

8 



s=0 



e— n 

t— 2 



«2-2s 



n 



+ — (1 -2s)£ 



-2..s 



s=0 



(3.10) 



The first sum is absolutely convergent, so one can simply set s = and evaluate it nu- 
merically. The second sum can be evaluated using zeta-function regularization. The result 
is 



- 2 E 

c ~ 2 



n- 



rr 



n(n — l)(n + 4) 
24 



(3.11) 



See Table [l] for evaluated for the first few lowest values of n. These results agree with 



those of |24|, where the sum (3.9) was regularized by a different method. 
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4 Beyond leading order: The strategy 

The 1/Nf corrections to the ground state energy on S 2 x R come from performing the 
functional integral over the gauge field. The effective action for the gauge field fluctuations 
obtained after integrating out the fermions is 



- I d 3 r d 3 r' y/g(?jA i {r)K n ' ii {r, r')Af{r') + ... 



(4.1) 



where the ellipsis denotes terms that are higher order in the gauge field fluctuations, and 
K n, ^(r,r') is the two-point correlator of the U(l) current J % : 



K n > ij (r,r') = -{J\r)Ji(r'))i 



(4.2) 



Here, i and j denote frame indices in the frame defined in (2.4). The subscript n on the 



angle brackets serves as a reminder that the expectation value should be evaluated in the 



monopole background (2.3). As can be read off from the action (2.1), the current J 1 is 



Nt 



(4.3) 



a=l 



In computing the two point function (4.2 ) we should treat the Nf fermions as non-interacting, 
so 



KUr,r') = -N f (^{r)a^{r)^{r')a^{r')) 



(4.4) 



where if> corresponds to a single charged fermion in the background gauge field given by 



(2.3). In terms of the fermion Green's function G n (r,r') = (if)(r)if>'(r')) n , one can write the 



kernel (4.4) as 



K%{r,r') = -N f tr {a t G n (r,r')a 3 Gi(r,r')) 



(4.5) 



The reason why one can ignore the higher order terms in (4.1) now emerges. The factor 



of Nf can be absorbed by rescaling the gauge field fluctuations Aj —> Ai/y/Nf. After this 
rescaling, the higher order terms in the gauge field fluctuations that were omitted from 



(4.1) become suppressed at large Nf, so the leading contribution to the S 2 xK free energy 
comes from the Gaussian integral over the gauge field fluctuations computed using just the 



9 



quadratic action (4.1). 



The action (4.1) can be written in diagonal form by expanding both the fluctuations 
Ai and the kernel in plane waves in the ffi. coordinate and vector spherical harmonics 
on S 2 . Just as in the case of the spinor harmonics discussed in the previous section, the 
vector harmonics can be defined by simultaneously diagonalizing the angular momentum 
operators S 2 , L 2 , J 2 , and J3. See Appendix [a] for more details. These vector harmonics 
are constructed from the usual S 2 spherical harmonics Y^ m because the gauge field does not 
experience monopole flux. They have spin s = 1, orbital angular momentum £, and total 
angular momentum quantum numbers given by 



(in 



V, 



Cm 



w, 



£m 



— 1 > , rrij 



3 = l 
j = £>l, rrij 



m . 



J 



+ 1 > 1 



m . 



m . 



(4.6) 



We thus expand the gauge field fluctuations as 
A\r) -- 



2^ 



00 i—i 00 1 

e=l m=-£+l 1=1 m=-l 

00 i+1 

e=o m=-l-l 



(4.7) 



where the summation ranges follow from (4.6). The kernel KL(r, r'), seen as a 3 x 3 matrix 



K 9 (r, r'), should also be expanded in terms of the vector harmonics as 



K n (r, r') 



00 1 



1 m=-e 



xK"M 



t fat a>\ 



where, if t > 0, K"(w) is a 3 x 3 hermitian matrix whose entries depend on £ and u. When 
£ = 0, Vi m and W^i) m do not exist, and this matrix should be thought of as 1 x 1. Note that 
the matrix K™(w) acts trivially on the frame indices of the harmonics and can be computed 



10 



by inverting (4.8 ): 



47T 



21 + 1 



<i 3 r \/ g(r)i 



E 



^(Vl)m( 



xK n (r, r') 



(4.9) 



r'=0 



where we made use of the symmetry under S 2 rotations and translations along R. Using 



(4.7) and (4.8), one can now write down the effective action (4.1) in almost diagonal form: 



dco x ^ x ^ 
2^^ ^ 

£=1 m=-< 



CO 



.Ira | 
l W 1 



Co) 



Am, 
ly I 



(4.10) 



The Gaussian integral over the gauge field fluctuations then gives, roughly. 



^^/^f>+l)logdetK^ 

•> 1=0 



CO 



(4.11) 



This expression is rough because it ignores a very important subtlety: gauge invariance. 
Indeed, the effective action Sa should be independent of the pure gauge modes that are part 



of the expansion (4.7), so it must be true that the matrix K"(u>) has many eigenvalues equal 
to zero. (To be precise, for each I there should be one eigenvalue that vanishes.) While it 
should be possible to work carefully in a fixed gauge, the subtleties related to gauge fixing 
are exactly the same for all n and they disappear from the difference^] 



i^-^i/^E^+Dlogdet 

1=1 



K?M 
K°( W ) 



(4.12) 



provided that we only take the ratio of the non-zero eigenvalues of these matrices. One can 
furthermore assume that Fq ^ = 0, because when n = the S 2 ground state energy equals 
the scaling dimension of the identity operator, which vanishes. 



It can be checked explicitly using the formulas in the following sections that K 



n.UVi 



[CO 



6 The sum starts at I = 1 because, as mentioned above, when £ — the matrix K™(w) has only one entry, 
which vanishes for all n because it corresponds to a pure gauge mode. 
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K e {uj) 



0, which implies that the matrix K"(u;) takes the form 

K 





K t (w) 




















{uj) \ 



K, 





n,WW 



(4.13) 



The entry K™ ,vv {uj) is a (non- vanishing) eigenvalue of this matrix. Since one other eigenvalue 
vanishes and since the trace of this matrix is the sum of all three eigenvalues, the third 
eigenvalue must be Kf uu + Kf . Let's denote the non-zero eigenvalues by 



j&* = KT VV 



(4.14) 



The labels E and B stand for ".E-modes" and "5-modes," respectively, because when re- 
stricted to S 2 one can check using (A. 10) that Ue m and Wt m have vanishing curl (just like 



an electric field), while Ve m has vanishing divergence (just like a magnetic field). 

To sum up, the correction to the free energy obtained after performing the Gaussian 
integral over the gauge field fluctuations is 



(n) 



doj 
2^ 



]T(2£+l)log 



K? E (u)K? B {u) 



uj) 



(4.15) 



Explicit formulas for K™' e {uj) and K™' b (uj) will be given in Section [fjj To find them, one 
should start by computing the fermion Green's function needed in (4.5), which is the goal 
of the following section. 



5 Fermion Green's function 



The Green's function in the background (2.3) is defined as 



G n (r,r') = (^W(r'))„. 



(5.1) 



12 



It satisfies (0 + J\) (r)G n (r,r') = -S(r — r'), so we can write it using the spectral decom- 
position as 



oo e-i 



-iu(t— t') 



m=— £ 



(5.2) 



where the matrices M n ^ and N Wj £ were introduced in (3.5). Using (3.4), it is not hard to see 



that acting with {ilj) + J(j (r) on the Green's function results in a delta function. 
To use this expression further, note that 



1 



nJ. 



lo + iM n/ co 2 + £ 2 - n 2 /A 
One can perform the u integral: 



N n /. 



00 ^ 1 ' -E n ,\T-r'\ 



ie 



S„ 



nlm \ 



X 



senfr — r')N 



-1 

1 



si, JO' A') 



where in order to simplify the subsequent expressions we defined the energy 



(5.3) 



(5.4) 



E n<e = v^ 2 -n 2 /4. 



(5.5) 



One can also perform the sum over m using the explicit expressions (A. 7) for the spinor 



harmonics. To understand the result, it is simpler to first consider a similar sum over m for 
the scalar monopole harmonics: 



) = e- m0 F. ,( 7 ) 



m=— I 



,'2e + l„ , s 2£ + l 

F& £(C0S7) = \ XTk £ _a(7, 0) = rr- 

>' 4^ 2>*> 2 v / 47r2 n / 2 



l + cos 7 ) ri/2 P^" ) (cos 7 ) 



(5.6) 
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where 7 is the angle between two points on S 2 , 



cos 7 = x ■ x' — cos 9 cos 9' + sin 9 sin 9' cos( 



(5.7) 



and the phase factor e 40 is 



e» oos^ = cos -cos - + 6-*"-" • 
2 2 2 



0' 
sin - sin — 
2 2 



(5. 



The addition formula (5.6) is the generalization of the notion of zonal spherical harmonics 



to the case of non- vanishing monopole flux. Note that when n > 0, the sum over m yields 
an expression that, up to a phase, depends only on the relative angle between the two points 
on 5* 2 . Fortunately, this phase will not play a role in the computation below. 



The addition formulas that generalize (5.6) to the case of the spinor harmonics are more 



complicated because of the spinor indices, and they will not be reproduced here. They give 



ic. 



-En,e\T-T'\ e -inB 



Q l nA c °si) (x-x')-a + sgn(r - r') 



x I Q^(cos7)l + Q^(cos7) (x + x) ■ a + iQ^(cos7)(x x x') ■ a 



(5.9) 



where the coefficients Q l n ^{x) (where x = cos 7) can all be expressed in terms of the function 



QnA X ) = 



1 



2£ + l 



F n Ax) 



1 



U- 1 



F n ^i(x) 



(5.10) 



as 



Qn 



E, 



Q 2 



x 



x 



n.C 



QnA x ) 



1 — X 

Ti 

— Qn,t{ X ) ) 



Ql 



IX 



Ql 



: Qn/ 

n 



x 



\x 



2(1 + a; 



QnA X ) 



Explicit care must be taken for £ = n/2 where one should replace (5.10) by 

(1 + xW 2 



Qn,n/2{x 



2 2+ ^7T 



(5-11) 



(5.12) 



In other words, only the first term in the expression for Q n /(x) in (5.10) should be considered 
in this case. The reader is referred to Appendix [B] for a check that when n = 0, the Green's 
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function (5.9) agrees with what one would expect from a conformal transformation of the 



flat space Green's function. 

Before moving on, let's discuss a few properties of the function Q n ,e{x) that will be 



important later on. In terms of Jacobi polynomials, the definition (5.10) reads 

(1 + x) n ' 2 



Q n ,e( x ) 



(0,n) 



(5.13) 



Using an identity that relates the difference of two Jacobi polynomials to a third Jacobi 
polynomial of a different rank, Q n ,e{ x ) can also be expressed as 



Qn,i( X ) 



\l+xT'\l-x) (1 , n) 



(5.14) 



Quite remarkably, the Jacobi polynomial in (5.14) can be traded for a scalar monopole 
harmonic with n — 1 units of gauge flux: 



Qn t e(cos9) = -> 



cos I 



4vr(£ 2 - n 2 /4) 



&Yn-X 



1 n+l 
" 2 ' 2 



(5.15) 



Because of this relation, the function Q n ,e{x) satisfies a second order differential equation 



that follows from the eigenvalue equation (|A.4|) for the monopole harmonics: 

QnA X ) 



2 1 + i 







(5.16) 



This equation also holds in the case £ = n/2 where one should use (5.12) instead of (5.10). 

Since the function Q n ^ arose from describing fermions with total angular momentum 
j — £ — 1/2, it may not be very surprising that it can be expressed in terms of a monopole 
harmonic with angular momentum equal to j. It is less clear why the scalar monopole 



harmonic that appears in (5.15) experiences one fewer units of magnetic flux than the spinor 



harmonics. In particular, when n — 1, Qi/{x) can be expressed in terms of the usual 
spherical harmonics Y+~ , or equivalently in terms of the associated Legendre polynomials 
PHx). 
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6 Evaluating the gauge field kernel 



6.1 General formulas 



Using (4.9)-(4.14) and (5.10)— (5.11 ), as well as the explicit formulas for the spinor and vector 

x i 



harmonics, one finds 



Kl' B (u) 



(4 



7T 



E n j> + E ni 



(2£ + !)£(£ + 1) u 2 + (E n4 , + E nr y 



(4tt) 



(2f+l)f(l + l)S W H(^l%) : 



E n j' + E ni 



2 



+ E n! ffE n; e"Z 2 



+ Ent'Ent'iIn 



(6.1; 



with 

x i 
X E 



i 

2 I dx 



2£(e + l)F 0>e + (l-x)F^ 
1 +x 



Q 2 Qn,e'Qn,£" — (1 — x2 )Eo,lQ'n,£>Q'r 



nl" 



(6.2) 



1 + x 

— j dx FQ zQn^iQnju , 

J -i v — x 

Zf = 2 y dx (Fq^ - (1 - x)F^) (g 2 Q n ^Q n ,,» - (1 + xfQ'^Q'^,,) , 

Xf = y dx (F^ + (1 + X)F^) QnfQnf' ■ 

Since both F ^ and Qn i obey second order differential equations and since all the expressions 



in (6.2) are symmetric under exchanging £' and £", (6.2) can be brought into the canonical 



form 



dx [a(x)F 0/ (x) + (3(x)Fq / (x)] Q n ,t<{x)Q n ^{x) 



(6.3) 



with some functions a(x) and (3{x) that depend on which expression in (6.2) we are consid- 
ering. This canonical form is useful when evaluating the integrals with respect to x. 
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6.2 Particular cases 



6.2.1 n = 

Let us start by evaluating the kernel K^' (u). By "inspection" of the integrals in (6.2) for 



many values of £, £', and £", it is possible to guess the general formula 

(f + £"){{£' - t"f - £ 2 ){{£' + £"f - £{£ + !)){£ + £' + £" - 1) 



—y 



2tt£(£ + 1) ff n (u 2 + (£' + £"f)(£ + £> + £") 

(i'-l £" -\ £-\\ 

x 



£" - 1 £ - 1 




(6.4) 



where factor on the second line is the square of a 3-j symbol. At fixed £', the summation over 
£" is of course convergent because the 3-j symbol vanishes when £" > £' + £ — 1. However, the 
remaining summation over £' is divergent because the terms in the sum approach the constant 
value — l/(47r) at large £'. This divergence can be regularized by adding and subtracting 
1/(47t) from each term and then using the zeta-function identity YH»=i 1 = C( — 1) — ~~ 1/2: 



00 / 1\°°1 1°°/ i \ 

E «»H + - E s = s + E + S J ■ (6.5) 

£'=1 v 7 £'=1 t>=\ v 7 



The summation over £' is now convergent and can be performed explicitly at fixed £. For 
instance, for £ = 1 we have 



Similarly, 



K»' E (n\ 1 ,V 4f 2 (l + u; 2 )-a, 2 ^cothf_ , fi fi , 

^ M = 8^ + 47r (4f 2 _ l)( 4 f 2 + u, 2 ) = (2 + ^ } 16(1 + u; 2 ) • (6 - 6) 



of, x / 9N (1 + w 2 )tanh^ 
A# M = (6 + w 2 ) V J 



n p 9x w(4 + w 2 ) coth ™ 



16 W (4 + ^) (67) 



2 
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and so on. It is not hard to see that in generaQ 

£(£+l) + u 2 



(6.* 



where 



D° e (co) 



T{{£ + l + iu)/2) 



AT((£ + 2 + iuj)/2) 



(6.9) 



is the scalar field kernel appearing in eq. (35) of |33|. 

One can perform a similar analysis for K 0,b (uj). By inspection of (6.1)-(6.2), I found 



T A),B, \ 



2ir£(£ + 



—y 



(£' - £") 2 (£' + £")(£ +!-£'- £")(£ + £' + £") 
u 2 + (£' + £") 2 

?' - 1 £" -I 



(6.10) 



x 







0. 



The sums over I' are again divergent but can be regularized precisely as in (6.5). A few 
particular cases give 



0)B (1 + W 2 )tanh^ 



16u 

tsO,b, s w(4 + w 2 )coth 

K 2 M = 



2 



K 3 (w) 



16(1 + UJ 2 ) 

;i + w 2 )(9 + w 2 )tanh 



(6.11) 



7TUJ 
2 



16cu(4 + lo 2 ) 

and so on. From these expressions one can guess the general formula 

£ 2 + io 2 



K„ (lo) 



(6.12) 



The quantity that appears in (4.15) is the product K®' e (lo)K®' b (lo), which is given by 

£(£+ 1) +u 2 



256 



(6.13) 



7 The same expressions for (w) and K e ' (w) were obtained in collaboration with Subir Sachdev as 



part of a similar computation in the CP^ model 
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6.2.2 n = l 



When n — 1, the expressions (6.2) written in the canonical form (6.3) are 



Xf = {Ji - J 2 



I2 — — J\ — J 2 , 



+ !)-£ 



12 pill 



+ 



+ 1) [£{£ + 1) J - 2 J 2 ] - [Ji - J 2 + + 1) J ] 



p/2 _j_ £//2 _ 



(6.14) 



Xf = J x + J 2 -£(£ + l)J 



where Jj are integrals involving Fqj, Qi^{x), and Qi/"(x): 



J {£,£',£" 



1 1 

F ^(x)Qi t£ '(x)Qi t £»(x) 

1 1 X 



7-1 1 — X 

J 2 {£, ?,£") = [ dx-^F^(x)Q 1>e ,(x)Q ir (x) 
J-i l — x 



(6.15) 



One needs to evaluate these integrals. In doing so it is helpful to recall that (5.15) implies 
that Qi/(x) is expressible in terms of the associated Legendre polynomials: 



QiA x ) 



V2£ 



Vl -xPl_ 1/2 {x) . 



n(2£- 1){2£+ 1) v ~ —«-i/2V 
We should also recall that Fqi is proportional to a Legendre polynomial: 



2£ 4- 1 
F 0/ (x) = ^±±P,(,) 



(6.16) 



47T 



(6.17) 



When substituting (6.17) into Ji and J 2 we should also make use of the following formulas 



relating the derivative of the Legendre polynomials to associated Legendre polynomials: 



P' l {x)= l -[Pl l {x)+t{t + l)Pl l {x 

xP^x) = 1 -[p?(x) + e(e + i)p°(x)} . 



(6.18) 
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The J; then become integrals over products of three associated Legendre polynomials; they 
can be evaluated using the formula 



dxP^(x)P^(x)P^(x 



-1 







4 4 4 
mi m 2 -m 3 , , , 



n 



(4 + Q! 
(ii-Trn)!.' 

(6.19) 



which holds for m 3 = mi + m 2 . The resulting expressions for J + , J_, and Jo are relatively 
messy, but each of these quantities can be expressed in terms of 3-j symbols that can be 
easily evaluated using a computer program: 



Jn 



1 

2 



,"+4,^+4 



1 

2 



(4 + |) (4 + |) (2£ + 1) (l 4 4 4 £ 2 



167^4(4 + 1)4(4 + 1) \0 



.0 1 



, (6.20) 



■h 



(4 + |) (£2 + \) (2£ + l)y/£(£ + 1) [ I + 1 4 4 





32ttV4(4 + 1)4(4 + 1) 



x 



y^ + 2)(£ + 3) 



+ 1 



£+1 



1 -1 
(6.21 



J 2 



{£i + \) (4+20(21 + 1)^1(1 + 1) (t 4 4 

32vrV4 (4 + 1)4(4 + 1) VO 



s/(e-i)(e + 2) 



I tl t 2 
-2 11 



+ 1 



1 



(6.22) 



With (6.20)-(6.22) and (6.14), one can now evaluate K\ ,e (uj) and K\' h '(u) using (6.1). I 



checked that plugging (6.20)-(6.22) into (6.14) yields the same answers as performing the 



integrals in (6.2) explicitly. 



In (6.1), the sums over £' and £" run from 1/2 to infinity. While there are no problems 

i" > 1/2, extra care must be taken when £' — 1/2 or £" = 1/2. 



with using (6.14) when 



It is not hard to see using (5.12) that the contribution to K\' E (u) from £' = 1/2 and/or 
£" — 1/2 vanishes, while the contribution to Kl ,B (u) is given by 



4ttuj 2 + £(£+1 



(6.23) 
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So one can restrict the sums in (6.1) to run over £',£" > 1/2 and add (6.23) to K^' B (uS). 



6.3 Numerics for n = 1 



One can calculate the corresponding contributions to the free energy (4.15) when n = 1: 



^/|X>+l)>o g § 



Kl' E (u) 



(6.24) 



These expressions are defined so that the total 0(N®) correction to the free energy is 
F{ + F{ . Imposing a relativistic cutoff on the summations and integrations in (6.24) 
by restricting I and u in the range 



+ 1)+lo 2 <L(L + 1) 



(6.25) 



for some cutoff energy scale L, I find that both F{ and F± diverge logarithmically with 
L: see Figure [l] However, as can be seen in Figure |2j the logarithmic divergence cancels out 

F (1)E p(l)B 

0.6 
0.4 
0.2 

1/L 



-0.2 
-0.4 
-0.6 
-0.8 



0.05 



0.10 



0.15 



0.20 



Figure 1: The quantities F^ E (brown) and F[ 1>a (red) evaluated with the relativistic cutoff 



(1)B 



(6.25). They each diverge logarithmically as L — >■ oo. 



from the sum F^\ 

Extrapolating to L = oo, I get 



(i) 



-0.0383. 



(6.26) 
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0.034 



0.035 



0.036 



0.037 



0.038 




0.05 



0.10 



0.15 



0.20 



1/L 



Figure 2: The correction = F[ 1)E + F[ l)B to the free energy evaluated by summing up 



the expressions in (6.24) with the relativistic cutoff (6.25). The orange points are obtained 



by evaluating (6.24) numerically, and the solid line is a cubic fit. 



This is the sought-after correction in (2.6) to the ground state energy on S in the presence 



of one unit of magnetic flux through the two-sphere. 



7 Discussion 



The main result of this paper is that the scaling dimension of the monopole operator that 
inserts one unit of magnetic flux in QED3 with Nf fermion flavors is 



0.265iV7 - 0.0383 + 0(1/JV» 



(7.1) 



which is obtained by combining the leading order result 24 given in Table [X] with ggjjfr . The 



0(1) correction in (7.1 ) was obtained by performing a Gaussian integral over the fluctuations 
of the gauge field on S 2 x R around a background constant magnetic flux of 2tt that is 
uniformly distributed throughout the S 2 . 

A striking feature of this computation is that the fermion Green's function on S 2 x R in 
the presence of n units of magnetic flux (as well as the kernel that appears in the quadratic 
action for the gauge field fluctuations, which can be written in terms of the fermion Green's 
function) can be expressed in terms of the spinless S 2 monopole spherical harmonics of [35p6] 
in the presence of n — 1 units of magnetic flux. This major simplification is what made the 
computation in this paper at all possible, and it would be interesting to understand its 
origin from a more conceptual point of view. In the case n = 1 analyzed here, the fermion 
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Green's function is particularly simple because it is expressible in terms of the usual spherical 
harmonics on S 2 , or equivalently in terms of certain associated Legendre polynomials. 

Another feature of the computation presented in this paper is an exact cancellation of 
UV divergences. Because of symmetry under rotations on S 2 and translations along R, 
the quadratic action for the gauge field fluctuations is diagonalized by a combination of 
vector spherical harmonics on S 2 and plane waves in the K direction. There are two gauge- 
invariant sectors depending on the properties of the vector spherical harmonics on S 2 : one 
sector involving vector harmonics with zero curl (the ".E-modes"), and one sector involving 
vector harmonics with zero divergence (the "5-modes"). The contribution of each of these 
two sectors to the S 2 ground state energy can be evaluated independently. Strikingly, each 
contribution is logarithmically UV divergent, but the divergences cancel exactly when the 
.E-modes and the 5-modes are added up together. This cancellation is far from obvious at 
intermediate stages of the computation, and it therefore provides a check on the method 
used here. 



A few comments on the significance of the result (7.1 ) are in order. It can be noticed that 



the 0(1) term in (7.1) is small compared to the coefficient of the leading term proportional 
to Nf. Since the expansion (7.1) is likely to be an asymptotic series, the smallness of the 
0(1) term suggests that the approximation (7.1) might be accurate even for small Nf. It is 
then reasonable to use this approximation to extract estimates for the upper limits on the 
number of fermion flavors below which compact QED 3 confines and below which the naive 
CFT limit e 2 — > oo in non-compact QED 3 might break down. 

As mentioned in the introduction, in compact QED3 monopole operators can proliferate 
and lead to confinement provided that they are relevant in the RG sense, i.e. if their scaling 



dimensions are smaller than 3. From (7.1 ) it is easy to see that the n = 1 monopole operator 
is relevant if Nf ^ 11.47 and irrelevant otherwise. I would therefore expect that compact 
QED3 should be in a deconfined phase precisely if Nf > 12. This result is consistent with the 
lattice analysis of 37 , which shows that the monopoles proliferate for Nf = 2 and Nf = 4, 
as well with the results of 16 , where it is argued based on the F-theorem 12 15 and the 



work of Vafa and Witten 38,39 that confinement would not be possible for Nf > 12. 

Lastly, one should stress that the computation presented in this paper was performed 
under the assumption that QED3 flows to an infrared conformal fixed point obtained by 



taking e 2 — > 00 in the action (1.1). This assumption is certainly correct at large Nf, where 
the whole RG flow can be studied perturbatively. As discussed in the introduction, below 
a critical value of fermion flavors, Nf < Nf nt , one expects the infrared physics to be signif- 
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icantly different from what the CFT limit e 2 — > oo would predict, the scenario supported 
by lattice data being that of spontaneous chiral symmetry breaking. For Nf < Nf Vlt it is 
reasonable to speculate that something would go wrong with taking the limit e 2 — > oo in 



(1.1). Prima facie evidence that this naive CFT limit is no longer appropriate would be 
that it leads to operator dimensions that violate the unitarity bound. Such a situation is 
not uncommon. In supersymmetric gauge theories, it is well-known that as one decreases 
the number of flavors it may happen that at some point certain operators hit the unitarity 
bound; to continue decreasing the number of flavors, one must pass to a dual description 
of the CFT because the original one is no longer valid. In non-supersymmetric QED 3 such 
a dual description may not be readily available, and for fewer fermion flavors new physics 
would be expected. 

One can therefore estimate Nj nt by assuming that for Nf < Nj vlt the naive CFT limit 
predicts unitarity bound violations. The unit strength monopole operator is a Lorentz scalar, 



so in a unitary theory its scaling dimension should be no smaller than 1/2. Using (7.1), one 
finds a unitarity bound violation for Nf < 2.03. Of the local operators that are polynomials 



in the fundamental fields, the one with lowest scaling dimension is ipl^a, for which 17 
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[*a] = 2-^- + 0(l/iV / 2 ). (7.2) 

The scaling dimension of this operator violates the unitarity bound if Nf < 1.44. While one 
should treat these bounds very cautiously, they suggest that the monopole operators are more 
constraining. That the monopole operators are more important at small Nf is consistent 



with the similar findings of 26,40 in supersymmetric theories. Lastly, it is intriguing that 
the estimate for Nf Tlt obtained this way is consistent with the lattice results, which suggest 
there is chiral symmetry breaking at Nf = 2, but not for larger values of Nf. 



Acknowledgments 

I am extremely grateful to Subir Sachdev for collaboration on related projects and for many 
insightful conversations. I also thank Liang Fu, Patrick Lee, Stefan Meinel, and Senthil 
Todadri for helpful discussions, as well as Ethan Dyer, Igor Klebanov, and Mark Mezei for 
their interest in my work and for their comments on preliminary versions of this paper. 
This work was supported in part by a Pappalardo Fellowship in Physics at MIT and in part 
by the U.S. Department of Energy under cooperative research agreement Contract Number 
DE-FG02-05ER41360. 



24 



A Monopole harmonics via conformal transformation 



A.l Scalar harmonics 

The (scalar) monopole spherical harmonics Y q 



of 35,36 can be used as a basis for 



the angular dependence of a charged scalar field in a background with n = 2q units of 
monopole flux sitting at the origin of M. 3 . Just like the usual spherical harmonics, they are 
eigenfunctions of the angular momentum operators L 2 and L3 with eigenvalues given by 



L 2 Y q/m = £(£ + l)Y qtim , L 3 Y q/m = mY q 



q,£m 



(a.i; 



If the electromagnetic potential is given by (2.3), the angular momentum operators can be 
written as 



L-a = -id^ - q , 

r 2 _ V72 , 2 <? 



sin 2 6 



(cosfl - 1)L 



(A.2) 



3 • 



where V 2 is the usual Laplacian on the two-sphere. It will be useful for us to write 



Y 



(A.3) 



The L 2 eigenvalue equation can then be written as 

-d x [(1 - x 2 )d x o q/m (x)} + m2+ f_ + J qmx Q q M*) = *(e + i)e,^(x) . (a.4) 

The regular solution of this equation can be given in terms of the Jacobi polynomials: 



! I (2£ + !)(£- m)\(£ + m)\ (1+x)^ ( _™- m) 



(A.5) 



The normalization factor in (A.5) is chosen such that \Y q ^ m (8 , <f))\ 2 integrates to unity over 
the two- sphere. 

Let's denote the scalar monopole harmonics Y q ^ m as \£, m), suppressing for now the extra 
label q. They transform in the spin-£ representation of the 5*0(3) rotation group. In these 
conventions, q is allowed to take half-integer values, and £ — \q\ > is an integer. 
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A. 2 Spinor and vector harmonics 

A basis of functions for a charged field of spin s in the same monopole background can 
be constructed using the angular momentum addition rules. The (2s + 1) components of 
this field can be written as \s,m s ), with — s < m s < s, and they transform in the spin-s 
representation of 5*0(3). A basis of monopole spinor harmonics is given by 



\£,s,j,mj) = (-1) 



-l+s— m 



s J 



x rrij — m s m s — m, 

m»=— s \ J J / 



\£,rrij - m s ) <8> \s,m s ) , 
(A.6) 



where I wrote down the Clebsch-Gordan coefficients explicitly in terms of the 3-j symbols. 
In particular, for s = 1/2 there are two sets of modes 



T, 



g,lm \ 



'22 2 



' e+m+l y 
2£+l I q/rn 



S, 



q,£m \ 



1 1 

m H — 
2' 2 



£—m 



y 

2£+l 1 qAm 



1 l+m+l y 

21+1 Y qAm+l){ a : 




(A.7) 



(This is what's called 



Vj±l/2,jm 



in (24).) Note that for T q ^ m we have — i — 1 < m < £, and 



for S q ^ m we have — i < m < i — 1, as follows from the fact that in both cases —j < rrij < j. 
If the monopole harmonics Y q ^ m are normalized such that they have unit norm, then T q ^ m 
and S q / m will also have unit norm. 

Similarly, when q = 0, one can define the vector harmonics corresponding to s = 1: 



u tm (e,<i>) = \e,i,£-i,m) , 
v tm {e,<f>) = \e,i,e,m) , 

W tm (9,<f>) = \£,lJ + hm) . 
In thinking about vector harmonics on S 2 it is customary to define 



1 



1 



-AY, 



tm \ 



(A.8) 



(A.9) 



where the indices i,j = 8, 0, and e 9 ^ = 1 is the unit antisymmetric tensor. A straightforward 
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calculation shows 



1 



Y lm {0,4>)dT + 



2£+l 



2£+l 



- Xi 



tm\ 



(A.10) 



£ + 1 
2£ + l 



B A check: Fermion Green's function at n = 



As a check, when n = 0, the Green's function should be 



G {r,r>) 



i \x\ \x'\ a ■ (x — x') 



47T \x — x' 



3 ' 



(B.l) 



where x was defined in (2.4) in terms of the coordinates on S x R. Apart from the \x\ \x'\ 



factor, this expression is just the fermion Green's function in flat space, and \x\ \x'\ is the 
conformal factor needed to map the theory from IR 3 to S 2 x R. We can rewrite Go as 



G n (r, r') 



i a • (e(— ')/ 2 x - e^'-^x') 



(B.2) 



4vr (2 cosh(r - r') - 2 cos 7) 3/2 ' 
Let's now try to reproduce this expression from the spectral decomposition used in Sec- 



tion 5 When n — 0, (5.4) becomes 



oo £-1 



=1 m=-l 



(B.3) 



The spinor addition are 



m=-£ 
€-1 



E So,tm(8, 0)^o, (£-1 



(T 



47T 



(B.4) 



[— xP £ '(cos7) + x'P^'_ 1 (cos7 
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Then 



G o(r, r') = ; a 



Air 



[-x'P;(cos 7 ) + xP^cosj)] (-e(r' - r)) 



+ [-£P;(cos 7 ) + x , P;_ 1 (cos 7 )] 0(r - r') 



(B.5) 



Next we should use the generating function for the Legendre polynomials: 



y/l - 2tx + t 2 



:b.6) 



n=0 



Differentiating with respect to x and setting x = cos 7 and t — e ' r r '', one obtains 

P \r-r'\/2 



(2 cosh(r — r') — 2 cos7) 3 ^ 2 , „ 



£p;(cos T )e-^' 



(B.7) 



Using this formula in (B.5), one further obtains 



G (r, r') 



47r (2 cosh(r — r') — 2 cos 7) 3 ^ 2 
•f ' " 2 £ + e 



(T 



- e \T-rV2 jjt , + e -\T-TV2 £]d(r ,_ T] 



x r — r 



(B.8) 



Combining the terms in the square brackets, the final result is 



Gn(r, r') 



z 



)/2£ _ g(r'-r)/2^/ 



4vr (2 cosh(r - r') - 2 cos 7) 



3/2 



(B.9) 



which agrees with (B.2). 
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